Voltage-induced Shapiro steps in a superconducting multi-terminal structure 
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When a superconducting tunnel junction at a finite voltage is irradiated with microwaves, the 
interplay between the alternating Josephson current and the ac field gives rise to steps in the dc 
current known as Shapiro steps. In this work we predict that in a mesoscopic structure connected to 
several superconducting terminals one can induce Shapiro-like steps in the absence of any external 
radiation simply by tuning the voltages of the leads. To illustrate this effect we make quantitative 
predictions for a three-terminal structure which comprises a diffusive superconductor-normal metal- 
superconductor (SNS) junction and a tunneling probe, a set-up which can be realized experimentally. 

PACS numbers: 74.45.+c,74.50.+r,73.23.-b 



I. INTRODUCTION 

In 1962 Josephson predicted that a tunnel junction 
between two superconductors at a finite bias V would 
sustain an alternating current with frequency 2eV/h, 1 
Josephson also argued that the alternating current would 
be frequency-modulated by an applied rf field and this 
would lead to steps in the dc current-voltage characteris- 
tics at bias voltages given by V = mhv/2e, where v is the 
frequency of the field and m is an integer. In 1963 Shapiro 
shone microwaves onto a junction and observed the pre- 
dicted steps^ which since then are referred to as Shapiro 
steps. This experiment constituted the first proof of ac 
Josephson effect and also paved the way for applications 
of this effect £ 

The Shapiro steps are a consequence of the interplay 
between an ac Josephson current and a microwave sig- 
nal. The observation of Shapiro steps does not necessar- 
ily require an external microwave generator: when two 
Josephson junctions are electromagnetically coupled like 
in the classical experiment of Giaever^ a voltage-biased 
junction can be used as a microwave source and induce 
Shapiro steps in the current- voltage characteristics (I-V) 
of the other one. More subtle coupling schemes have 
been achieved, in which two superconducting weak links 
placed a few micrometers apart share a common elec- 
trode. Shapiro-like steps were then observed in the TVs, 
which were attributed to the phase-locking of the Joseph- 
son frequencies when the dc voltages on two bridges 
are approximately matched. 5 Although there is no rig- 
orous microscopic theory of this non-equilibrium effect, 
it is usually attributed to the diffusion of the quasiparti- 
cle charge imbalance generated in and around the weak 
links £ Very recently, it has been predicted that one could 
also generate Shapiro steps by coupling a tunnel Joseph- 
son junction to a mechanical oscillator. 7 

In this work we predict that a different type of Shapiro- 
like steps can appear in a single multi-terminal supercon- 



FIG. 1: (Color online) Schematic representation of the setup: 
a metallic wire (N) is connected at its ends to superconduct- 
ing reservoirs Si and S r , biased at potentials — U/2 and U/2, 
respectively. An additional superconducting probe (S p ) at a 
potential V is attached to the wire through a tunnel junction. 



ducting structure without any external radiation. The 
idea goes as follows. Inspired by the experiment of 
Ref. [8j], let us consider the three-terminal structure 
shown in Fig. [TJ It is formed by a diffusive SNS junc- 
tion, and a superconducting tunneling probe S p attached 
to the normal wire. If one applies a potential difference U 
across the wire, ac currents with a frequency 2eU/h and 
its harmonics will flow along the N wire. If a potential 
V is now applied to the probe electrode, it will generate 
ac currents with Josephson frequencies 2e(V ± U/2)/h. 
The interference between these ac currents which both 
involve the quasiparticles in the diffusive wire gives rise 
to contributions to the dc current through the interface 
N-Sp, which in this setup appear as Shapiro steps at dis- 
crete voltages V — mU/2, where m is an odd integer. 
Notice that in this case the steps appear due to the non- 
equilibrium state created in the diffusive wire. In this re- 
gion the energy distribution function of quasiparticles is 
driven to oscillate at different frequencies by the dc volt- 
ages of the leads and the Shapiro steps are a result of the 
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beating of these internal Josephson frequencies. As illus- 
trated below, this effect is described quantitatively using 
the quasiclassical theory of superconductivity for diffu- 
sive systems, in particular time-dependent Usadel equa- 
tions. Finally, it is important to emphasize that these 
voltage-induced Shapiro steps may appear in a great va- 
riety of mesoscopic structures and, as explained below, it 
constitutes a valuable way to directly test the ac Joseph- 
son effect. 

The rest of the paper is organized as follows. In section 
II we remind the basic physics of diffusive SNS junctions 
at finite bias. Section III is devoted to the description 
of the theoretical formalism and contains the main tech- 
nical results. In section IV we discuss the numerical re- 
sults for both the current-voltage characteristics of the 
tunneling probe and the dependence of the height of the 
voltage-induced Shapiro steps on the different system pa- 
rameters. Finally, in section V we briefly summarize the 
main conclusions of this work. 



II. BASIC PHYSICS OF DIFFUSIVE SNS 
JUNCTIONS 

The effect described above is a consequence of the non- 
equilibrium properties of diffusive SNS junctions. The 
physics of these junctions is the result of the interplay 
between the proximity effect in the normal wire and the 
occurrence of multiple Andreev reflections (MARs). The 
proximity effect is the modification of the properties of a 
normal metal in contact with a superconductor and it has 
been extensively studied in diffusive hybrid nanostruc- 
tures. 9 MARs are the tunneling processes that dominate 
the subgap transport in SNS junctions. Here, successive 
Andreev reflections at both S electrodes lead to a pro- 
gressive rise of the quasiparticle energy until its energy 
exceeds the gap energy, Ad The interplay between prox- 
imity effect and MARs in diffusive SNS junctions gives 
rise to a rich variety of phenomena. Thus for instance, the 
conductance exhibits a very peculiar subgap structured 
Dubos et alH studied the ac Josephson effect by shin- 
ing microwaves onto a junction, and they observed the 
appearance of fractional Shapiro steps, which is a sig- 
nature of a non-sinusoidal supercurrent-phase relation. 
On the other hand, using the setup of Fig. [T] Pierre et 
al& measured the non-equilibrium distribution function 
in a long silver wire, where the proximity effect was neg- 
ligible. They showed that this function exhibits several 
steps, which is a manifestation of MARs. 

The consequences of the interplay between proximity 
effect and coherent MARs in the dissipative current has 
been addressed theoretically only very recently by us- 
ing the time-dependent Usadel equations^ This theory 
describes the crossover from the short junction regime 
(L £ = \J hD / A),— where L is the wire length, £ is 
the superconducting coherence length and D the diffu- 
sion constant, to the incoherent limit (L ^> £)d The 
intermediate regime is the relevant one for the observa- 



tion of the voltage-induced Shapiro steps. 

III. CALCULATION OF THE 
VOLTAGE-INDUCED SHAPIRO STEPS 

We consider the structure depicted in Fig. [IJ where 
the position of the tunneling probe along the wire is de- 
noted by x (xi = and x r = L). We assume that the 
three superconductors (S;, S r and S p ) have the same en- 
ergy gap A. Our goal is the calculation of the current 
I{V) through the probe junction S p when a voltage U 
is applied across the wire. We assume the N-S; jr inter- 
faces to be fully transparent and neglect phase-breaking 
phenomena. 

To solve our problem we use the quasiclassical theory of 
superconductivity for diffusive systems.— This theory is 
formulated in terms of momentum averaged Green func- 
tions G(R, t,t') which depend on position R and two 
time arguments. These propagators are 2x2 matrices in 
Keldysh space ("), where each entry is a 2 x 2 matrix in 
electron- hole space ("): 

*■-(££)■ <» 

The Green functions for the left (1) and right (r) leads and 
for the probe (p) electrode can be written as Gj(t,t') = 
e -t^(i)f 3 /2fiQ o ( t _ tl y^(t')H/2h^ where ig the 

phase of the order parameter of the electrode j = l,r,p 
given by d<pj(t)/dt — 2fij/h, where the chemical poten- 
tials are \n = —eU/2, [i r — eU/2 and \x v — eV. Here, 
Go(t) is the equilibrium bulk Green function of a BCS 
superconductor. We now transform to energy representa- 
tion, in which the propagator G(R, e, e') depends on two 
energy arguments. It satisfies the non-stationary Usadel 
equation, which in the N region reads 

(G o VG) + ef 3 G - Gf 3 e' = 0, (2) 

where T3 is the Pauli matrix in electron-hole space. The 
convolution product o is defined as (A o B)(e,e') = 
f dei A(e, ei)B(ei, e'). Equation ([2]) is supplemented by 
the normalization condition 6 o G = — 7r 2 <5(e — e')I. Due 
to the finite bias U, to solve Eq. ([2]) is a formidable task, 
which is explained in detail in Ref. What matters 

for our discussion is that the Green functions in the wire 
adopt the form: 

G(R, e, e') = ]T G , m {e)5(e m - e'), (3) 
in 

where e m = e + meU. Here, m is an even integer 
for the diagonal components of the Green functions in 
Nambu space and an odd integer for the off-diagonal 
ones. Eq. © means, in particular, that the energy dis- 
tribution function and the density of states in the normal 
wire oscillates with the Josepshon frequency 2eU/h and 
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its harmonics. These oscillations result in a parametric 
pumping of the N-S p junction and hence to Shapiro steps. 

Assuming that the electrode S p is weakly coupled to 
the N wire, we can express the time-dependent current 
I(V, t) through the tunnel probe up to first order in the 
tunneling conductance Gt as follows^ 



I{V,t) 



= (-) 

\87re J 



J cftiTr{f 3 [G tu (t,t 1 ),G p (t 1 ,t)] A '} 



(4) 

where G w is the Green function of the wire at the posi- 
tion of the probe junction and G p is the Green function 
of the probe electrode, i.e. a bulk BCS Green funtion. 
Using Eqs. (|3l4p it is easy to show that 



w*)= E E WK (n 



<p+2neVt/H+meUt/h) 



(5) 



m— — 00 n— — 1 



where </> the dc part of the phase difference. From Eq. ([5]), 
one can distinguish two contributions to the dc current, 
which we shall simply denote as /. First, there is a back- 
ground current Ig = Iq that contributes at every voltage 
V . More important, there is a series of contributions at 
discrete voltages V m = mU/2 with m odd^ which are 
given by 7 sha p iro (0) = E^oi 1 !™^ + I^e'^iV - 
V m ), where = Thus, Shapiro-like resonances 

are induced every time the probe voltage V is an odd mul- 
tiple of U/2. From Eq. ([4}, one can show that the ac 
current components in Eq. ((5|) can be expressed as 



1 l — 



( G 7 



r de {[^) , m (^p^m + eV) 

J — OO 

+ [T«] Qm (e)^(e m + eV)}.(6) 



Here, (e) are Fourier components of the 

anomalous Green functions induced in the normal wire 
by proximity effect. These components contain the infor- 
mation of ac Josephson effect in the SNS junction and de- 
pend on L, U and x. Obviously, when L>( these com- 
ponents vanish and in turn the Shapiro stepsJ^ The func- 
tions Tp' K {e) are the bulk Green functions of the elec- 
trode S p given by F£{e) = -ttA/^A 2 - (e - i0+) 2 and 
pK(e) = -2ilm{^ p 4 (e)}tanh(/3e/2), where (3 = l/k B T 
is the inverse of the temperature. 

For completeness let us also mention that the back- 
ground current can be approximated by 



hi = ( ~ ] / < !< <"/' ( ' 



eV) Pw (e) [f( e -eV)-f w (e)} 



where /(e) is the Fermi function, f w (e) is the dc part of 
the distribution function of the N wire at the position of 
S p , pp(e) is the bulk BCS density of states and p w (e, U) — 
Im{ \Gw\ 00 }/it is the non-equilibrium spectral density in 
the normal wire. 




FIG. 2: (Color online) Zero-temperature dc current in the 
probe electrode S p at zero bias (U — 0) for different lengths. 
The probe is located on x = L/2. The inset shows the critical 
current as a function of the wire length. 



IV. RESULTS AND DISCUSSION 

It is illustrative to start the discussion of the results 
with the case of zero-bias (U = 0). In this limit there are 
no Shapiro steps, but the dc current exhibits a supercur- 
rent peak at V=0. This can be seen in Fig. [21 where we 
plot the dc current through a probe located on x = L/2 
for different wire lengths^ Apart from the supercurrent, 
the main feature of these I-V curves is the presence of a 
gap, which is equal to A + A 5 , where A g is the well- 
known minigap present in the N wire due to proximity 
effect (see for instance Fig. 1 in Ref. 13). Such a mini- 
gap decays with L as A g ~ 3.2eT, where = hD/L- is 
the Thouless energy. In the inset of Fig. [2] we show how 
the critical current decays with L. This critical current 
is a good reference for the height of the Shapiro steps at 
finite U. Notice in particular, that for L <C £ the critical 
current adopts the value I = 7rGTA/2e, which is sim- 
ply the critical current of a tunnel junction between BCS 
superconductors i^i 

Now, we turn to the case of finite bias U. To give a first 
impression of how the steps would appear in an experi- 
ment, we show in Fig. (21a) the total contribution to the 
dc current, i.e. background current plus Shapiro steps, 
for L = 2£, x = L/2 and three different values of the volt- 
age U. Notice that since we are assuming that the system 
is voltage biased, the Shapiro steps appear as peaks in 
the current, rather than steps as in a current-biased con- 
tact.— We have also plotted in Fig. [3] both the distribu- 
tion function and the non-equilibrium spectral function, 
which determine the shape of the background current. 
The first two peaks in the I-V at V = ±U/2 could sim- 
ply be viewed as the supercurrent peaks that correspond 
to the condition of equality between the potentials of the 
probe Sp and of one of the reservoirs S; !r . However, their 
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FIG. 4: (Color online) Height of the Shapiro steps at zero 
temperature as a function of the position of the probe (x) for 
eU = 0.5A and different wire lengths. 



FIG. 3: (Color online) (a) Zero-temperature dc current in 
the probe S p located on x = L/2 for L = 2£ and three dif- 
ferent values of the voltage U: from bottom to top eU/A = 
0.3, 0.6, 0.9. (b) dc part of the distribution function of the nor- 
mal wire at x = L/2 as a function of the energy for the three 
voltages above, (c) Corresponding non-equilibrium spectral 
density in the normal wire (x = L/2) as a function of the 
energy. In the three panels the two upper curves have been 
shifted upwards for clarity. 



height depends on U and they are smaller than the crit- 
ical current at U = because the quasiparticles in the 
wire are out of equilibrium. Notice, also, the large sub- 
gap current, which is absent in Fig. [2] More obvious 
manifestations of the interference of ac Josephson cur- 
rents are the higher-order Shapiro steps, which appear 
at voltages that do not correspond to the alignment of 
the potential of S p to that of a reservoir. In Fig. [3ja), 
they are clearly seen at V = ±3t//2 at the lowest volt- 
age (U — 0.3 A), but they progressively disappear as U 
increases. 

Let us now study systematically the best conditions 
for the observation of the voltage-induced steps. First, 
we analyze their dependence on the wire length and the 
position of the probe at a given voltage U. In Fig. 0] we 
show for eU = 0.5A the height of the steps, denoted as 
S m u/2 for a step at V = mU/2, as a function of x for 
different wire lengths. Notice that we show the steps for 
positive voltages V and the first one for negative voltages, 
S—U/2- Indeed, one can show that in this setup the fol- 
lowing relation holds: S- m u/2{x/L") — S m jj/2(1 — x/L) 
for m > 0, which can be seen in the upper panels of 
Fig. [U The most remarkable features of Fig. [J] are the 
following. First, the step Sjj/2 is of the same order as 
the critical current at U = 0, although a bit smaller. 
Second, due to the biasing in this setup (see Fig. [TJ, the 
step 5(7/2 vanishes close the left electrode. The same hap- 



pens with the steps S3U/2, SbU/2, ■■■ at both boundaries 
because we are assuming that the electrodes are ideal 
bulk reservoirs. Third, the higher-order steps (m > 1) 
are, for this particular voltage U, at least an order of 
magnitude smaller than SV/2- Finally, with respect to 
the length dependence, the steps progressively disappear 
as L increases. For a given position, they decay roughly 
as £/L for L 3> £, which is the usual decay of induced su- 
perconducting correlations in a normal diffusive system. 
However, for L of the order of £, the dependence is not 
neccesarily monotonous, as it can be seen for the steps 
S3U/2 and S 5U / 2 - 

The other crucial parameter that controls the height of 
the steps is the voltage U. As Eq. © shows, in order to 
have a step at V = mil /2 (m odd), one needs to have a 
non-zero component ] 0m - Such an ac component 

is related to the coherent transfer of n = (\m\ + l)/2 
Cooper pairs through the SNS junction, which requires 
the occurrence of coherent MARs of at least order n. 
An rt-order process contributes significantly for voltages 
2A/n < eU < 2A/(n— 1) and its probability decreases as 
the voltage increases. Thus, one naively expects (i) high- 
order (m > 1) steps to be more clearly visible at low bias 
U and that (ii) at voltages U > 2A, only the steps S±u/2 
survive. Indeed, these expectations are confirmed by the 
calculations, as can be seen in Fig. [5l where we show 
the height of the first two steps (for positive voltages) 
as a function of U for different wire lengths and a probe 
located on x = L/2. Notice that the step Su/2 reaches 
its maximum at ~ 2A and then decays very slowly as 
the voltage increases whereas S 3 u/ 2 vanishes rapidly for 
voltages above the gap. 
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eU/A eU/A 

FIG. 5: (Color online) Height of the first two Shapiro steps 
for positive voltages in x = L/2 at zero temperature as a 
function of U for different lengths L. 



induced steps are a direct manifestation of the ac Joseph- 
son effect in the SNS junction and they are visible as long 
as there are sizable superconducting correlations induced 
in the N wire. We stress that these steps can be induced 
in any coherent mesoscopic structure attached to several 
(more than two) superconducting electrodes. Moreover, 
if in the setup of Fig. [1] all the electrodes were strongly 
coupled to the central wire, new phenomena would oc- 
cur like the appearance of fractional Shapiro steps for 
V = (m/n)U/2, with m,n integers. Thus, the effect pre- 
dicted here suggests the appearance of new physical phc- 
nonema and in particular, provides a new way to directly 
test the ac Josepshon effect in mesoscopic structures. 



V. CONCLUSIONS 

In summary, we have predicted the possibility to gen- 
erate parametrically Shapiro steps in multi-terminal su- 
perconducting structures by tuning the voltages in the 
reservoirs. We have illustrated this effect with detailed 
calculations for the case a diffusive SNS system with a 
third electrode coupled to the N region through a tunnel 
junction. Our results, based on the solution of the time- 
dependent Usadel equations, show that these voltage- 
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